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Comment on "Density Matrix Renormalization 
Group Study of the Haldane Phase in Random 
One-Dimensional Antiferromagnets" In a recent 
Letter Hida presented numerical results indicating 
that the Haldane phase of the Heisenberg antiferromag- 
netic (HAF) spin-1 chain is stable against bond ran- 
domness, for box distributions of the bond strength, 
even when the box distribution stretches to zero bond 
strength. The author thus concluded that the Haldane 
phase is stable against bond randomness for any distri- 
bution of the bond strength, no matter how broad. This 
conclusion contradicts our earlier prediction that a 
sufficiently broad bond distribution drives the HAF spin- 
1 chain through a critical point into the Random Singlet 
phase. In this Comment, we (i) point out that the ran- 
domness distributions studied in Ref. ^ do not represent 
the broadest possible distributions, and therefore these 
numerical results do not lead to the conclusion that the 
Haldane phase is stable against any randomness; and (ii) 
provide a semiquantitative estimate of the critical ran- 
domness beyond which the Haldane phase yields to the 
Random Singlet phase, in a specific class of random dis- 
tribution functions for the bond strength. 

(i) The box distribution function for bond strength 
J used in Ref. [Q takes the form Pbox{J) = for 
1 - WI2 < J <l + W/2, and Pi,ox{J) = otherwise. 
There is no randomness for = 0, and the random- 
ness increases with W . The maximum randomness for 
this class of distribution functions is reached a.t W — 2, 
as further increasing W would introduce ferromagnetic 
bonds to the system and change the physics completely. 
On the other hand, the W = 2 distribution does not 
represent the broadest possible distribution among all 
possible bond distributions. In fact, it is a member of 
the class of power-law distributions: Ppow{J) = aJ~^^" 
for < J < 1, and Ppow{J) = otherwise. The case 
= 2 of Phox corresponds to a = 1 in Ppow The 
power-law distributions are of particular interest and im- 
portance in studies of random spin chains, because un- 
der real space renormalization group (RSRG) transfor- 
mations, all generic distribution functions will flow to 
power-law distributions in a HAF spin- 1/2 chain [^,1). 
The broadness of a power-law distribution is parameter- 
ized by the single positive exponent a > 0. The smaller 
a is, the broader the distribution is on a logarithmic scale 
(the width is I /a) Clearly, for any < a < 1, the 
distribution function Ppow{ J) is broader than those stud- 
ied in Ref. [|l| . Thus the numerical results of Ref. |^ do 
not necessarily lead to the conclusion that the Haldane 
phase in HAF spin-1 chain is stable against any bond 
randomness, no matter how broad the distribution Q. 

Hida also argues that the persistence of the string order 
is a consequence of the absence of translation symmetry 
in the effective spin-1/2 description of the random HAF 
spin-1 chain similar to what happens in the random 
dimerized spin-1/2 chain 0. We point out that the Ran- 


dom Singlet phase in the effective spin-1/2 model does 
not have translational symmetry, and that broken trans- 
lational symmetry and string order are not synonymous. 

(ii) While a naive extension of the original form of the 
RSRG cannot be applied directly to the random 

HAF spin-1 chain unless the initial distribution function 
is sufficiently broad, we have developed a RG scheme 
that is suitable for the random HAF spin-1 and other 
random spin chains. The key ingredient of this scheme is 
to project out the highest energy state(s) in the strongest 
bond and keep all lower energy states at each RG step. 
Details of this scheme will be presented elsewhere . Us- 
ing this scheme to simulate the HAF spin-1 chain with 
power law distributions, we have found two phases sep- 
arated by a critical point ac ~ 0.67. For a > ac (which 
includes the flat distributions of Ref. ||l|) the system is 
in the Haldane phase characterized by the string order, 
while for a < ac it is in the Random Singlet phase. While 
there are approximations involved and our estimate of ac 
may not be very accurate numerically, it clearly suggests 
that one needs bond distributions that are broader than 
those studied ^ to destabilize the Haldane phase. 
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